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Measurements were made on refractive index changes with hydrostatic pressures 
between 1 bar and 1 kbar using the helium yellow line. The materials studies were: KBr 
NaCl, LiF, diamond, MgO, quartz, AI2O3, and three silicate glasses. All the materials 
increased in refractive index with pressure except diamond, MgO, and A1 2 3 which decreased, 
and LiF which did not change. The results were compared with photoelastic measurements, 
and PockeFs geometric theory of photoelasticity was substantiated as well as Mueller's 
physical theory. The data show that the ratio of change of polarizability with density is 
greater for solids having stronger interatomic repulsive forces. Volume and temperature 
coefficients of polarizability were evaluated for the cubic crystals and glasses. The thermo- 
optic behavior of crystals and glasses is discussed. 



1. Introduction 

The relationship between index of refraction and 
density is of importance because of the effect of 
atomic interactions on the atomic polarizabilities. 
For transparent solids, values of pdp/dp, where p is 
the index and p the density, have been obtained by 
calculation from the photoelastic constants [l]. 1 For 
crystals as well as glasses values of pdp/dp so cal- 
culated usually have been found to be less than the 
corresponding values calculated on theoretical 
grounds, i.e., from the Lorentz-Lorenz and the 
Drude relationships [1]. The most extreme example 
of this discrepancy is probably found in the case of 
diamond and MgO where pdp/dp values calculated 
from photoelastic constants are negative, a result 
which is clearly impossible to reconcile with the 
assumption of constant polarizability involved in 
the initial derivation. The theory of photoelasticity 
due to Mueller [2, 3] recognizes the change in the 
intrinsic polarizability of an atom due to strain. 

For glasses, data on pdp/dp have been obtained by 
Ritland [4] (and others) in the annealing tempera- 
ture range where both index and density change 
with time at constant temperature and on volu- 
metric relaxation following removal of high hy- 
drostatic pressures [5]. In the former case pdp/dp is 
found to be somewhat less than theory predicts and 
in the latter case somewhat greater. 

Direct experimental data on pdp/dp are required 
to understand thermo-optic properties. Values of 
dp/dT, where T is the temperature, have been 
measured for many solids and it is known that the 
sign may be either positive or negative. For most 
crystals dp/dT is found to be negative but in the case 
of diamond, MgO, and ZnS the reverse is true [6]. 
Recent work has shown that for A1 2 3 dp/dT is 
positive for both the ordinary and extraordinary 



1 Figures in brackets indicate the literature references at the end of this paper. 



rays [7]. Most optical glasses exhibit a minimum 
in the p-T curve somewhat below room temperature 
so that at and above room temperature dp/dT is 
positive but at some lower temperature the sign 
reverses [8]. Any attempt to explain such a complex 
thermo-optic behavior must of necessity require an 
understanding of the pdp/dp effect because thermal 
dilatation is an important contributing factor to 
the value of dp/dT. 

Finally, by reversal of the customary procedure, 
data on pdp/dp may be used to evaluate photoelastic 
constants in certain instances. For glasses only two 
constants are required, so that knowledge of pdp/dp, 
plus information on the optical path difference pro- 
duced by a unidirectional stress serves to determine 
completely the photoelastic constants. 

For these several reasons it was considered of inter- 
est to obtain direct measurements of the density 
coefficient of refractive index, pdp/dp, in solids at 
constant temperature by the application of hydro- 
static pressure. This is the first time that such an 
approach has been used. This report contains data 
on pdp/dp obtained in the pressure interval between 
1 bar and approximately 1000 bars for several crys- 
tals and glasses of interest. Data were obtained 
from specimens of the cubic ionic crystals, LiF, KBr, 
NaCl, and MgO; the covalent cubic crystal, dia- 
mond; and the birefracting crystals, A1 2 3 and quartz 
where both the ordinary and extraordinary rays were 
studied. The glasses studied were fused silica, a 
commercial plate glass and a borosilicate crown glass, 
BSC 517/645. Measurements reported here were 
made only at 25 °C using the helium yellow line. 

2. Experimental Method 

In these experiments interference fringes were ob- 
served in a plate of the material having plane, pol- 
ished, nearly parallel faces [9]. The fringes were 
viewed in reflection using; collimated helium light of 
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\= 5875.62 A at normal incidence. For a specimen 
of thickness, t, and index, /x, the fringe number N is 
given by the expression 



N\=2tn. 



(1) 



Application of hydrostatic pressure will produce a 
change in /x, a decrease in t, and a shift in the inter- 
ference fringes. By differentiating (1) and rearrang- 
ing it is found that 



AfJL = 



AN\ At 



(2) 



A measurement of the number of fringes that pass a 
point of reference, and the thickness together with 
data on the initial index and the linear compression 
At/t permits calculation of the change in index. 

The experimental apparatus and method have been 
described earlier in detail and will only be outlined 
here [10]. Briefly, a plate of the test specimen is 
immersed in liquid in a pressure vessel equipped with 
glass windows. Hydrostatic pressure is generated by 
compressing the liquid and the number of interference 
fringes passing a reference line, i.e., the fringe count, 
was determined visually. Estimates of the number 
of fringes were made to the nearest tenth of a fringe. 

The change in fringe number observed in approxi- 
mately 1 kbar is not very large for these solids 
because of their rather low compressibilities com- 
pared to the liquids previously studied and because 
the decrease in thickness under pressure acts in such 
a direction as to decrease the fringe number (see 
eq (2)). In the most favorable case AN is known 
only to three significant figures. Extreme precision 
in the other quantities is not required. The most 
reliable values for the initial indices were used, 
however, and values for the index at the frequency 
of the He yellow line were generally obtained by 
interpolation of the dispersion data reported. The 
data on the indices of refraction of the solids may 
be found in references [11] through [20]. For BSC 
517/645 the index was reported for the sodium D 
line only. However, from the low dispersion of 
such glasses and the small frequency separation of 
the sodium and helium yellow lines, the value re- 
ported was assumed to apply to the present measure- 
ments within the required precision. 

Compression data for all materials except BSC 
517/645 have been reported, usually in terms of 
volume compression (AV/V ). Volume compres- 
sions were converted to linear values for the isotropic 

materials through the relationship A£/£o=7Tt7* For 

quartz and A1 2 3 , where data on linear compressions 
are available, At /to was computed for the direction 
perpendicular to the optic axis. The data for deter- 
mining the compression of the solids may be found 
in references [21] through [26]. (For BSC 517/645 
the linear compression was calculated from elastic 
constant data measured at 1 bar [27].) It has been 
estimated that, at most, the uncertainty in values of 
At 1 1 lies in the third significant figure. 



For the optically isotropic solids (cubic crystals 
and glasses) the measurements were straightforward 
but for the crystals of lower symmetry, A1 2 3 and 
quartz, which are optically anisotropic the method 
was modified to permit measurements for both the 
ordinary and extraordinary ray. For these materials 
the specimen was oriented so that the optic axis was 
perpendicular to the direction of the light beam. 
In the measurements a large polarizer was inserted 
between the source and the window of the pressure 
vessel. On rotating the polarizer two distinct fringe 
systems could be seen depending on whether the 
electric vector of the plane polarized light trans- 
mitted by the polarizer was parallel or perpendicular 
to the optic axis of the specimen. (For the quartz 
crystal which is optically active, it should be noted 
that, since the fringes were obtained by reflection, 
the emergent polarized light was still plane polarized 
in the same plane at which it was incident [28]. 
This would not have been true if the fringes were 
observed in transmission [16].) On applying hydro- 
static pressure the fringe shift for each separate 
fringe system was noted, giving data on both rays 
for the birefringent crystals. 

3. Preparation of Specimens 

Most specimens were prepared by sectioning or 
clearing larger samples of commercially available 
synthetic single crystals. However, both quartz and 
diamond were prepared from naturally occurring 
single crystals. The glasses were of commercial 
origin. Although the test specimen required was 
approximately 1 cm X 1 cm X 0.5 cm, it was found 
expedient to grind and polish larger specimens and 
then to cut the final specimen from the most perfectly 
ground portion of the large piece. Rough blanks of 
both A1 2 3 and quartz were cut witb the optic axes 
in approximately the correct orientation. The blanks 
were then oriented using the Laue back-reflection 
x-ray technique and ground so that the optic axis 
was in the surface of the specimen. 

3.1. Grinding and Polishing 

A brief account is given here of the preparation of 
the crystals because these techniques are not widely 
known. The specimens were first ground flat and 
with nearly parallel surfaces using a surface grinder 
for the harder materials and using hand grinding 
with fine abrasive for the alkali halides. The harder 
crystals were prepolished on a wood lap with 8^-22^ 
diamond powder using olive oil lubricant and then 
given a final polish on cloth impregnated with 3ju 
diamond dust using an alcohol-water solution as a 
lubricant. The alkali halides were given only the 
latter of the two polishing steps. A1 2 3 specimens 
were available only in small diameters and the 
grinding and polishing operation was conducted on 
an assemblage of five small disks, set in wax on a 
flat glass plate. 

The diamond was a 2 carat brilliant cut gemstone 
and the table face was found to be satisfactory for 
one surface. The culet was ground down parallel to 
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the table face using a high speed steel lap impreg- 
nated with diamond powder to produce a surface 
approximately 1 mm in width. Satisfactory fringes 
could be observed through this small surface. 

In all polishing operations it was necessary to 
check the fringe system at frequent intervals. If it 
was noted that the surfaces were deviating too much 
from parallelism adjustments were made in the 
polishing to reduce the angle between the surfaces. 

3.2. Aluminizing 

In order to improve the sharpness of the inter- 
ference fringes the reflecting power of the surfaces 



of the specimens was increased by vacuum deposition 
of a layer of aluminum. The front surface of each 
specimen was covered with a partly reflecting film 
while the rear surface was made fully reflecting. 
Multiple reflection fringes were obtained in this 
manner but because of the large thickness to length 
ratio of the specimens the fringes were not extremely 
sharp. A series of depositions of different thick- 
nesses were made on the front surface of plate glass 
specimens until the fringes appeared to reach a 
maximum intensity and sharpness. This film thick- 
ness was found to be satisfactory for all other 
specimens. 



Table 1. Pressure induced changes in refractive index as measured with heliu?n yellow radiation of 5875.62 A 



Material 


Refractive 
index, 

MHe 


Specimen 

thickness, 

k 


Pressure, 
P 


Thickness 
change, 
-A^oX103 


Observed 

change in 

fringe number 

-AJV 


Change in 

refractive 

index 

Am X 105 


KBr 


L.560 

1. 540 
1.392 
1.738 

2. 418 

1.545 
1.554 

1.769 
1.760 
1. 458 
1.517 
1. 518 


Cm 
0. 3210 
. 6748 
. 5683 
.3970 
. 4037 

.5112 
.5112 

. 5659 
.5659 
. 4967 
. 6693 
. 6566 


bars 
966. 1 
963. 6 
963. 2 

960. 8 
960.0 

957. 6 
961.5 

961. 8 
957.6 
960.2 
960. 1 
960. 


2.118 
1.348 
. 349 
.190 

.058 

.944 
.948 

.105 
. 105 
.865 
.717 
.723 


11.3 
21. 6 
9.4 
6.8 
3.4 

6.9 

7.5 

6.2 
6.2 
7.3 
8.9 
9.1 


227 


NaCl 


115 


LiF _. 





MgO___. 


-17 


Diamond- 


-11 


Quartz: 


103 




107 


Sapphire: 


-14 




-14 


Fused Si0 2 --. 


83 


BSC 517/645—. 


70 


Plate glass. .-_ _ 


69 







4. Results 



The results of the experiments are given in table 
1 which shows the decrease in fringe number, the 
dimensional change calculated at the pressure given 
and the change in index of refraction. It should 
be f noted that the change in fringe number is pre- 
sented in table 1 for the direction of increasing 
pressure, i.e., between 1 bar and the pressure listed. 
For birefracting crystals the index for the ordinary 
and extraordinary rays are noted by co and e respec- 
tively. Each value of AN given in table 1 represents 
the average of at least two determinations. Such 
duplicate measurements agreed to within ±0.1 
interference fringe. From the agreement between 
the duplicate measurements the limit of reproduci- 
bility in the index change is estimated to be 
±2X10~ 5 . 

From the tabular data it is seen that the change 
in index is quite small in absolute value, varying 
from in LiF to 2 in the third decimal place for KBr. 
The index decreases for MgO, diamond and both the 
ordinary and extraordinary rays of A1 2 3 . All 
other materials show an increase in index except 
LiF for which the index is unchanged. For the 
birefracting crystals both the ordinary and extraor- 
dinary rays show the same change in index within 
the limits of error. 



5. Discussion 

5.1. Cubic Crystals 

From the experimental values of Am and the com- 
pressibility data, it was possible to calculate values 

for the quantity p -~ where p is the density. For the 
cubic crystals studied here the theory of photo- 
elasticity shows that p —- can be calculated by the 
relationship 

(3) 



4ri (Pii+2p ^ 



where the elasto-optic constants, p ijf relate the state 
of strain existing in a crystal to changes in the index 
ellipsoid [6]. These values have been calculated for 
a number of crystals [1, 6] and the results are given 
in table 2 together with the initial indices, /x. [The 
values of these experiments obtained over the 
pressure interval of approximately 1 kb are denoted 
pA/*/Ap while those from the photoelastic theory are 
designated as pd^/dp]. It may be noted that the 
photoelastic data were obtained for the sodium D 
line and the present data for the helium yellow line. 
However, because of the small frequency difference 
between the two lines and the known low-frequency 
dependence of pdn/dp [29, 30] any differences arising 
from dispersion effects can be taken to be negligible. 
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Table 2. Values °f P\T~ ) an d A o for cubic crystals and 
glasses 





Refrac- 
tive 
index, 

Mo 


'(£) 


<$) 


L-L equation 


Drude equation 


Material 


>(t) 


Xo 


>(t) 


Xo 


KBr 

NaCl 


1.559 
1.544 
1.392 
1.736 
2.417 
1.458 
1.517 
1.518 


0.35 
.28 
.00 
-.31 
-1.58 
.32 
.33 
.32 


0.35 
.24 
.1 

-.40 

-.28 

.29 


0.68 
.66 
.44 
.97 

2.62 
.53 
.62 
.62 


0.48 
.58 
1.00 
1.32 
1.60 
.40 
.47 
.49 


0.46 
.45 
.34 
.52 

1.00 
.39 
.43 
.43 


0.24 
.38 


LiF.._ 


1.00 


MgO_ 


1.60 


Diamond 

Fused S1O2 

BSC 51 7/645—. 


2.58 
.17 
.24 


Plate glass. __ 




.26 









A comparison of the figures in table 2 shows that 
there is essential agreement between the predictions 
of photoelastic theory developed by Pockels [31] and 
the actual measured values. Considering the errors 
inherent in the measurements of the photoelastic 
constants [6], i.e., small fringe shifts, nonuniformity 
of stress distribution in uniaxial compression, etc., 
it is considered that the agreement in table 2 is 
satisfactory except in the case of diamond. For this 
material, however, the errors in the photoelastic 
experiments and in the present work are likely to be 
particularly large because of the small specimens 
available and the small changes that are observed. 
There is no question as to the negative sign of 
pA/x/Ap for diamond, and it is interesting to note that 
the results of the present experiments bear out the 
theoretical predictions of negative values of pA^/Ap 
for both MgO and diamond. 

The data for diamond obtained in this study are 
in much better agreement with the results of Gibbs 
and Hill [32] on the change of the dielectric constant 
of diamond with pressure. Gibbs and Hill noted 
that the dielectric constant of diamond, €=5.66 at 
27 °C at frequencies as high as 10 MHz is equal to 
the square of the optical index of refraction at infinite 
wavelength, juoo, reported by Peter [14]. Since p is 
proportional to 1/V, (d In P /dP) T = - (d In V/Z>P)t=P, 
where P is the pressure and /3 the compressibility. 
It follows that 



/dMo 



If 



2 

" Moo 5 



\dpj r \dP d 
it follows that 



dP 



P_\ = l/d/Ov 2 



1 /fri,\ 1 /be\ 
H m \bP/r 2e\dPj T ' 



(4) 



(5) 



If the value of p ( — J in the present study is used for 

p(d/x 00 /dp) T in eq (4), and employing the equivalence 
shown in eq (5), (l/e)(de/o\P) T for diamond has been 
found to have the value — 2.28X10" 7 bar _1 . The 
value obtained experimentally by Gibbs and Hill for 
a type I diamond is — 2.40X10" 7 bar -1 . 



2 The derivation leading to eq (4) was suggested by D. D. Wagman of NBS. 



5.2. Glasses 

Experimental values for pA/x/Ap for the glasses are 
given in table 2. The values are very similar for 
all three glasses and probably reflect the fact that 
the data are determined largely by the Si0 4 tetra- 
hedra common to all three glasses. Numerical 
comparisons with photoelastic theory can be made 
only for vitreous Si0 2 for which Jog [29] has reported 
the necessary photoelastic data. As seen in the table 
the value of pd^/dp calculated from the photoelastic 
constants is in reasonably good agreement with the 
value obtained here by direct measurement. All 
three glasses show positive values for pAju/Ap. 
Furthermore the photoelastic data on numerous 
optical glasses reported by Schaeffer and Nassenstein 
[33], Vedam [34], and Mueller [3] invariably indicate 
positive values for pdfx/dp. This fact will be discussed 
later. 

5.3. Noncubic Crystals 

It is interesting to note that the increase in index 
for both rays in quartz is greater than observed in 
fused silica. The result can be explained only by 
the conclusion that interatomic interactions play an 
effective role in the index change with density and 
that those are different in the crystal and the glass 
because of the crystal structure and the higher initial 
density of the crystal. From the data of table 1 it 
appears that the extraordinary ray suffers the larger 
change in index but the difference is probably within 
the experimental error. Pockels [35] has measured 
the photoelastic constants of quartz and his results 
can be extrapolated linearly to the present pressure 
of 957 bars to permit calculation of A/jl. Under these 
circumstances the photoelastic data yield calculated 
values for the change of index of 100X10 -5 for the 
ordinary ray and 104 X10 -5 for the extraordinary 
ray. These values are in good agreement with the 
results of the present studies. In contrast to quartz 
both rays of sapphire show the same change in index 
but both are negative. There do not appear to be 
enough photoelastic data available to permit calcu- 
lations of the index change for A1 2 3 [7). 

5.4 Photoelastic Constants of Glass 

From the present data on A/jl/AP and auxiliary 
data on the amount of double refraction introduced 
into glass by uniaxial compression it is possible to 
compute the constants q n and q l2 , where the piezo- 
optic constants, q ijy relate the state of stress which 
exists in a solid to changes in the index ellipsoid [6]. 
The relationships used [34] are as follows: 



and 



A M /AP=M 3 [<Zii+2fe]/2 
(A/i 2 -AMi)/AP=/x 3 [g 12 -2 11 ]/2, 



(6) 
(7) 



where, in addition to the terms previously defined, 
A/x 2 and A/xi represent changes in index for light 
polarized respectively perpendicular and parallel to 
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the stress direction. (In eq (7), AP denotes a change 
in a compression applied in one direction only.) 
The quantity (A/io — Ap.i)/AP is defined as the relative 
stress optical coefficient and has been reported to be 
2,81 brewstersforBSC 517/645 [36] and 2.62 brewsters 
for plate glass [37]. These data phis values of 
Aju/AP serve to determine values for r/ n and q V2 from 
eq (6) and eq (7). These values are given in table 3. 
Values of Neumann's photoelastic constants, p and 
q, are also given in table 3; the relationships used in 
evaluating them are given by Vedam [34]. 

The elasto-optic constants, p n and p u can also be 
calculated from the relationships given by Vedam 
[34]. Values of the elastic moduli needed in the 
calculations are given by Spinner for BSC 517/645 
[27]. From measurements made on commercial 
plate glass [38], values for Young's modulus, the 
rigidity modulus and Poisson's ratio have been 
found \o be 723.9 kbar, 297.9 kbar, and 0.215 re- 
spectively. These values were used to calculate the 
elasto-optic constants for plate glass shown in table 3. 

Table 3. The photoelastic constants of BSC 517/645 and 
plate glass 



Glass 


V 


Q 


Pu 


7>12 


f/iiX10' ;! 


tfl 2 X1013 


BSC 517/645... 


0.168 
.154 


0.087 
.087 


0.115 
.114 


0. 221 
.203 


cm 2 ldyne 

0.315 

.369 


cm 2 /di/ne 
1.92 


Plate glass 


1.87 



5.5. Polarizability and Density 

It seems self-evident that density changes must be 
accompanied by some change in molecular or atomic 
polarizability and it is useful to analyze the present 
data to ascertain what information can be obtained 
in this connection. The most comprehensive theory 
appears to be due to Mueller [2, 3] who considers 
the following factors involved in change of refractive 
index of a solid under applied stress: (1) change of 
density; (2) change of the Lorentz-Lorenz field; (3) 
change of the coulomb field; and (4) change in the 
intrinsic polarizabilities of the atoms. Mueller 
appears to be the first investigator to recognize the 
importance of the last factor. The theory has been 
applied so far only to cubic crystals and isotropic 
solids because in these cases the symmetry is such 
that the coulomb field can be taken to be zero. 
Under hydrostatic pressure the symmetry is un- 
changed and the coulomb field remains zero but the 
Lorentz-Lorenz field may change because of the 
change in density. 

There are two principal theoretical relationships 
relating the index and density, the Lorentz-Lorenz, 



and the Drude, 



M 2 -l/l\ 

7+iyp) =ai 

(m 2 — l)/p=«2 



(8) 
(9) 



where a x and a 2 are constants. In deriving the 



Lorentz-Lorenz equation the local field lias been 
accounted for through the theoretical value 4/7r3 P 
where P is the polarization of the medium and the 
change in local field with density is accounted for 
in the equation. The change in the local field with 
density has been ignored in the Drude equation. 3 

If a x and a 2 are rigorously constants eq (8) and 
eq (9) may be differentiated to obtain the change of 
index with density. This was apparently first done 
by Pockels [39] who obtained the two expressions 



p^=(m 2 -D(m 2 +2)/6 m 



5m 
dp" 



: IM 



-1)/2m 



(10) 



(11) 



corresponding to the Lorentz-Lorenz and Drude 
formulations respectively. In eq (10) and eq (11) 
partial derivatives are used to indicate that there 
may be a change in atomic polarizability which is 
neglected. Mueller [2, 3] introduced the strain 
polarizability parameter, X , to evaluate the change 
in polarizability by means^of the relationship, 



pdfi/dp= (1 — X )pd/x/dp. 



(12) 



This relationship follows from considering a, and a 2 
to be dependent on the polarizability. By using 
eq (12) in conjunction with eqs (10) and (11) and 
the experimental values of pA/i/Ap values of X have 
been calculated for both the Lorentz-Lorenz and the 
Drude formulations. These data are shown in 
table 2. 

Burstein and Smith [1] have proposed that X f) is 
a measure of the degree of homopolar binding. 
According to Burstein and Smith X () should attain a 
maximum value for intermediate bonding, i.e., 
bonding containing appreciable amounts of both 
ionic and homopolar bonding, and decrease as the 
bonding becomes either more ionic or more homo- 
polar. The present studies do not agree with this 
conclusion, the value of X being largest for the most 
covalent crystal-diamond. It appears that the 
conclusions of Burstein and Smith were based upon 
erroneous data for diamond. 

For the glasses, calculated values of X are also 
given in table 2. These values show little variation 
for the three glasses but it is noteworthy that X 
data calculated by both the Lorentz-Lorenz and the 
Drude equations are positive. Previous data on 
liquids show that X is positive when computed from 
the L-L equation but negative when computed from 
the Drude equation [30]. This behavior in liquids 
was explained as arising from the change in the local 
field brought about by compression of the liquids 



3 Krishan and Roy [40] have shown the equivalence of the Lorentz-Lorenz 
and Drude equations in expressing the frequency dependence of the refractive 
index and have shown that the Drude expression' accounts for the local field by 
incorporating its effect into the value of the fundamental al (sorption frequency. 
However, eq (9) is a simplified version of the generalized Drude equation and is 
valid only at a single frequency. When written in the form of eq (9) it is implied 
that the change in the local field with density has been ignored. 
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and incorporated into a change in the fundamental 
absorption frequency in the Drude model. This 
greater change in the local field of the liquids arises 
from their greater compressibilities. 

Ramachandran and Radhakrishnan [41, 6] have 
shown that the thermal properties of the index of 
refraction may be understood by considering the 
change of index with temperature to consist of 
three independent contributions which they label 
P, Q, and R. P represents the change arising only 
from a change in the number of scattering centers 
and is a pure density effect, Q represents the change in 
polarizability resulting from the density change, and 
R consists of the effect of temperature alone on the 
polarizability. Recognizing that dp/dT=P+Q+R f 
the values they derived using the Drude equation 
are as follows: 

P-v(^) 03) 

«~<»aHr) <"> 



R- 



(I/JL . dfJL 

~-dT +yp dp 



(15) 



where 7= 



Values of pd^jdp they used 



1 (<k\ 

p \dTJp 

[41, 6] were evaluated from the elas to-optic constants 
(eq (3)). The quantities P, Q, and R have been re- 



evaluated using the same data that Krishnan and 
Roy used except that experimental values of pA/x/Ap 
have been substituted for pdp/dp. [There is one 
further change; it appears that the earlier results 
were calculated using an erroneous value for the 
volume expansivity of diamond, 1.3X10~ 5 /°C [41, 6] 
compared to the correct value, 0.31X10 _5 / o C [42, 9].] 
The results are given in table 4. The present data 
agree reasonably well with the earlier results in 
most instances, there being slight changes in Q and 
R resulting from the present values of A/x/Ap. There 
is a large change, however, in the results for diamond. 
It is also possible to calculate R, which is defined 
as (dn/dT) Vj from constant volume data. For a 
given material the calculations involve finding a 
temperature at 1 atm where the specific volume is 
the same as that at 25 °C and an elevated pressure; 
the refractive indices which pertain for the two 
different situations may then be used to find R. 
Such calculations have been made using values of y 
and (dn/dT)p as a function of temperature. The 
original sources for the values of y and (dfx/dT) P for 
the different solids have been identified in the 
list of references. 4 Much of the necessary infor- 
mation may be found in the compilation made 
by Krishnan [6]. The values are given in the column 
headed R (experimental) in table 4. It can be seen 
that there is reasonably good agreement between 
corresponding values of R in the two columns. 



* BSC 517/645 is known to have the same nominal composition as the glass 
designated as BSC-2 in reference [8]. 



Table 4. Analysis of change in refractive index and polarizability for cubic crystals and glasses 



Material 


"d 


7X105 


dnldTXlW 


PX1W 


QX10* 


RXW 


#X105 
(experi- 
mental) 


/Vdct\ 
\ctdV/ T 


\adV Jt 


\aOV ) T 

X105 


(I *2) 

\adTjv 

X105 


\<xdTjv 

X105 
(experi- 
mental) 


KBr 

NaCl 

LiF 


1. 559 
1.544 
1.392 
1.736 
2.417 
1.458 
1.517 
1.518 


°c-i 

12.0 
12.0 
10.2 
3.9 
0.31 
.15 
2.0 
1.9 


°c-i 

-3.6 
-3.8 
-1.3 

1.6 
1.0 
0.92 
.17 
.27 


-5.5 

-5.4 

-3.4 

-2.3 

-0.3 

-.06 

-.9 

-.8 


1.3 

2.0 

3.4 

3.5 

0.8 

.01 

.2 

.2 


°C-i 

0.6 
-0.4 
-1.3 

0.4 
.5 
.97 
.8 
.9 


°<7-i 
0.4 
-.3 
-1.3 
0.6 
.6 
1.0 
.7 


0.2 
.4 
1.0 
1.5 
1.7 
0.2 
.2 
.3 


-0.8 
-.7 
0.0 
.5 
.7 
-.8 
-.8 
-.7 


OQ-1 

-9.4 

-7.8 

0.0 

2.0 

0.2 

-.1 

-1.5 

-1.4 


1.6 

-3! 8 

0.7 
.8 
2.5 
1.9 
2.0 


°C-i 

0.9 

-.7 
-3.9 


MgO 


1.0 


Diamond 

Fused Si0 2 

BSC 517/645 

Plate glass . . 


0.6 
2.6 
1.5 









For the present purposes where volume and 
temperature may both be considered as independent 
variables it appears more desirable to change the 
analysis slightly as follows: Starting with the 
Drude relationship (eq (9)) and considering the 
polarizability a 2 to be a function of V and T, it 
can be shown that 



dpjp,= 



(m 2 -i) 

2 M 2 



(16) 



(In eq (16) and the treatment that follows a 2 is 
written simply as a.) In this expression - ( ^ J is 



the pure temperature coefficient [of polarizability, 
— I ^ryr) is the pure volume coefficient of polariza- 
bility and the term ( — ( r-=, ) — 1 ) represents the 
total contribution of volume to the change in index. 
At present the term — ( yj/) can be evaluated from 
the constant temperature experiments. To obtain 

_ ( ^tft ) the volume is taken to be a function of J P 
a \Ol/v 



and T. Then dV= 



j dT-\-l Yp ) dP and using 
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U %1 



V\dTj P ~' and ~V\dPj T = 

cast in the form 



= (3 eq (16) can be re- 



dp 



M — 



2m ; 



l \}M\ dT+{ydT ~^ P) 

\a \dVj T 



'V/daX 



0) 



(17) 



The value of V/aCda/dV)^ determined from eq (16) 
can be combined with available data on y and on the 
temperature variation of index at constant pres- 
sure to derive values for - ( ^7^ ) • If the Lorentz- 

a \ol/v 

Lorenz equation is used the same expression will 
be obtained except that the factor - — 



replaces 



2m 2 



6m 2 
This represents only a different 



scaling factor and in no way affects the relative 
magnitudes of the terms in the brackets which are 
of primary interest. 

The calculated values of the terms of interest in 
eq (17) are given in table 4. The volume coefficient 
of polariz ability is seen to be positive in all materials 
as expected. The positive sign means a decreased 
polariz ability on compression and an increased 
polariz ability on heating arising from the thermal 
expansion. In the alkali halides the numerical 
values of the coefficient increase as the anion packing- 
increases. Values for MgO and diamond are the 
largest in the table, a result consistent with the large 
interatomic interactions presumed to occur in these 
materials. The small values for the glasses are to 
be expected from the smaller interatomic inter- 
actions expected for such open structures. The 
smallest value in the table applies to fused Si0 2 
which has the most "open" structure. (The figures 
in table 4 have been rounded off and this conclusion 
is not obvious from the data given.) The total 
volume contribution is given in the next column and 
these values are generally negative except for dia- 
mond and MgO. It is noteworthy that in every 
case the volume coefficient of polariz ability repre- 
sents an appreciable portion of the total volume 
effect. It must be concluded that any relationship 
between index and density which ignores the effect 
of this contribution to the index can, at best, be only 
approximately correct. The next column tabulates 
the total volume term multiplied by the expansivity 
to permit comparison on an equal basis with values 
for the temperature coefficients of polariz ability 
given in the next to last column. The temperature 
coefficients vary in both sign and size with the value 
apparently decreasing with packing in the alkali 
halides and being rather large positively in the 
glasses. There are insufficient data, however, to 

verify that these conclusions are typical. (~>^) 

may also be found by substituting values of (dfji/dT) v 



in eq (16). These values are shown in the last 
column of table 4. Again it may be seen that there 
is fairly good agreement between corresponding 
values in the two columns of figures. 

With the data given in table 4 it is instructive to 
consider the change in index resulting from two 
separate experiments, first an isothermal compression 
and second an isobaric temperature change. An 
isothermal compression produces an increase in index 
in KBr, NaCl, and all the glasses because the pre- 
dominant effect is an increase in the number of 
scattering centers per unit volume. The decrease 
in polarizability with decreasing volume is less than 
the density effect in these materials. In LiF the 
index undergoes no change (to 1 kbar) because the 
density increase and the polarizability decrease 
exactly balance. In both MgO and diamond the 
index decreases on compression, the polarizability 
decrease overbalancing the density increase. 

It is of interest to consider the effect of pressures 
higher than the 1 kbar used here. It appears likely 
that the volume coefficient of polarizability will 
increase as the volume decreases because the inter- 
atomic interactions should increase in magnitude. 
In the absence of other effects, this requires that at 
some elevated pressure the indices of the solids which 
increased with pressure at low pressures will pass 
through a maximum value and then decrease at still 
higher pressures. Further the decrease in index 
should accelerate at still higher pressures. The 
pressures required to produce these changes are not 
obvious but plans are being made to extend these 
measurements to 10 kb. In this range it appears 
very probable that a decrease might be expected 
for LiF. For the other alkali halides a reversal or a 
decrease in the rate of change of index with density 
might be expected. 

An isobaric temperature increase results in a 
decrease in indices for KBr, NaCl, and LiF but an 
increase in indices for MgO, diamond and all the 
glasses. Although the resultant effects are similar 
for members in each of the two groups, there are 
various causes, and it is necessary to consider the 
details of the process. In KBr the volume expansion 
causes a decrease in index because the number of 
scattering centers per unit volume decreases. This 
effect is great enough to overbalance the increased 
polarizabilities arising from the increased volume and 
temperature. In NaCl the temperature coefficient 
acts in conjunction with the density effect to decrease 
the index. In LiF there is no total volume contri- 
bution and the decrease in index arises solely from 
the temperature coefficient of polarizability. In both 
MgO and diamond the index increases because the 
positive temperature and volume coefficients of 
polarizability overbalance the negative contribution 
from the decrease in the number of scatterers with 
the volume coefficient playing the major role. All 
the glasses show an increase in index because the 
comparatively large temperature coefficient of polar- 
izability outweighs the negative effect of the total 
contribution of change in volume which is nearly as 
large. 
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The glasses require further consideration because it 
is known [8] that the index-temperature relationships 
of most optical glasses exhibit a minimum in index. 
The glass BSC-5 17/645 has such a minimum near 
— 80 °C and although data on plate glass do not 
appear to be available a minimum may reasonably 
be expected at some temperature below 25 °C. On 
the other hand the index of refraction of fused Si0 2 
is known to decrease monotonically with decreasing 
temperature [51,16] despite the reversal of its co- 
efficient of thermal expansion. A reasonably sound 
explanation of these effects requires data on the two 

terms - ^ T7 ^^ , and - ^7™ which are apparently not 
a Vol a 01 

available at this time but certain conclusions may be 
of interest until such data are obtained. It might 
be expected that for normal materials the polariza- 
bility arising from pure thermal agitation should de- 
crease at higher temperatures because of increased 
interaction. However, it is not apparent whether 

^-j, increases or decreases. Data by Radhakrishnan 

on LiF [49] show that d\x\dT becomes increasingly 
negative at higher temperatures. At room tempera- 
ture this crystal shows a zero volume contribution, 
and, if it is assumed that this contribution is negligi- 
bly small at elevated temperatures, it may be con- 
cluded that - ( rr-=£ ) is negative. The glasses, how- 
ever, are not normal solids so that it is not certain 
that they exhibit a similar behavior. It is known that 

7" ( VT* ) * s P os iti ye for the glasses and it is to be ex- 



V, 



V 



pected that — ( r—^ j would increase somewhat at 

low temperatures because of the closer packing. 

In table 4 it is observed that for BSC 517/645 and 
plate glass that the total volume contributions and 
the temperature coefficients of polarizability are of 
opposite signs but nearly equal in absolute values. 
If either one or both of those quantities changes 
slightly in the proper direction on lowering the 
temperature it is clear that the sign of dn/dT will 
reverse. Since this reversal occurs below room 
temperature in BSC 517/645 and is expected in 
a similar temperature range for plate glass, it is 
probable that the temperature coefficient of polariza- 
bility decreases with decreasing temperature in these 
glasses in contrast with the opposite behavior ob- 
served in LiF. The situation in fused Si0 2 is some- 
what different in that - ( ^-^ ) is at least ten times 
a \ol/v 

as great as the total volume effect although they are 
of opposite sign. In this situation a small change 
in either quantity will be insufficient to reverse the 
sign of d/jt/dT. In addition when the expansivity of 
fused Si0 2 becomes negative both the volume term 

and - -^-Tp are positive and cannot cancel. It may be 

noted that the interpretation here is similar to that 



advanced by Prod'homme [52] for the anomalous 
temperature-index behavior in optical glasses but is 
somewhat more detailed. It appears that the posi- 
tive value of djji/dT in optical glasses arises from the 
comparatively large temperature coefficient of polar- 
izability coupled with the low coefficient of thermal 
expansion. Since these quantities are not common 
to all glasses [52] it appears that they must be attrib- 
uted to the silica network. 
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